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Could Leptons, Quarks or both be Highly Relativistic Bound States of Minimally 

Interacting Fermion and Scalar? 
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The possibility that leptons, quarks or both might be highly relativistic bound states of a spin-0 
and spin- 1/2 constituent bound by minimal electrodynamics is discussed. Typically, strongly bound 
solutions of the Bethe-Salpeter equation exist only when the coupling constant is on the order of or 
greater than unity. For the bound-state system discussed here, there exist two classes of boundary 
conditions that could yield strongly bound solutions with coupling constants on the order of the 
electromagnetic fine structure constant a. In both classes the bound state must have spin one half, 
thus providing a possible explanation for the absence of higher-spin leptons and quarks. 



I. MOTIVATION FOR COMPOSITE MODELS 
OF LEPTONS AND QUARKS 

The motivation for studying composite models of 
leptons and quarks is almost entirely circumstantial. 
Three times in the past 150 years particles that were 
- or are currently - thought to be elementary could be 
organized into families: (1) In the 1860's Meyer and 
Mendelev arranged the atoms into the periodic table 
that consists of 18 families. The existence of fami- 
lies was ultimately explained by the realization that 
atoms are composite. (2) Almost 100 years later two 
of the then-fundamental particles, the neutron and 
proton - as well as many other baryons and mesons - 
were placed into SU(3) multiplets with the ultimate 
result that the neutron and proton were no longer 
viewed as being fundamental but instead were com- 
prised of quarks [E Hi- (3) By the mid 1970's many 
physicists realized that the charged leptons, the neu- 
trinos, the positively charged quarks and the nega- 
tively charged quarks constitute four families. Today 
the three charged families are known to have three 
members each, and the neutrino family almost cer- 
tainly has three members as well Q. Because the 
existence of families of "fundamental" particles has 
twice been explained by the realization that the par- 
ticles are actually composite, although speculative, as 
is any physics beyond the standard model, the most 
conservative approach to explain the existence of the 
lepton and quark families likely is to assume that the 
particles are composite. In the 1970's the hypotheti- 
cal constituents of quarks and leptons were given the 
name "preons" 

The circumstantial evidence for composite leptons 
and quarks is certainly not conclusively corroborated 
by direct experimental evidence. The factor relates 
the muon's magnetic moment fi and spin s according 
to the relation, 
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g^ = 2{l + a^). (1) 



The present experimental world average determina- 
tion for 15] is al^'P = 116592080(63) x 10~" . The 



primary source of error in the theoretical calcula- 
tion of flp arises from the contribution of the strong 
force, which cannot be calculated from first princi- 
ples but, instead, must be determined from exper- 
imental data. A direct determination of the con- 
tribution from strong interactions uses experimental 
data from electron-positron collisions, and the theo- 
retical value calculated from the standard model is Q 
^SM direct ^ (116591785 ± 51) x 10"". The theoret- 
ical value in is smaller than the experimental value by 
3.6 a. The indirect method [3| uses experimental data 
from hadronic tau decays, conservation of the vector 
current and isospin corrections and yields the value 
^SM indirect ^ (116591932 ±52) x 10"", which is 1.8 
a below the experimental value. The discrepancy sug- 
gests that there is physics beyond the standard model. 
But in addition to muon substructure, the discrepancy 
could also be explained, for example, by the existence 
of supersymmetric particles or by composite W and 
Z bosons. 

If leptons, quarks or both are composite, several 
characteristics that any composite model must possess 
can immediately be determined from the mass- and 
spin-spectra of the leptons and quarks. The masses of 
the charged leptons and quarks satisfy the inequalities 

nioloctron ^ mmuon ^ HIt- , mup ^ Ulcharm ^ nitop, 

and mdown < nistrange < mbottom- Siucc no structure 
has been conclusively detected for leptons or quarks, 
any composite system must be very strongly bound. 
This precludes a composite model whereby, for ex- 
ample, the electron, muon, and tau are successively 
much more massive because one or more of their con- 
stituents are successively much more massive. The 
large differences in the masses of the particles in each 
family are consistent with strong binding only if the 
bound system is highly relativistic. 

Even though the existence of families is explained 
by the atomic model of atoms, the quark model for 
mesons and baryons and, as discussed here, perhaps 
a preon model for leptons and quarks, each of the 
composite models is radically different. In the atomic 
model, aside from small mass defects, the mass of an 
atom is the sum of the masses of the constituents. For 
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mesons and baryons, most of the mass results from the 
kinetic energy of the quarks so the mass of a meson 
or baryon is substantiaUy greater than the sum of the 
masses of the constituent quarks. If leptons, quarks or 
both are highly relativistic bound states, the mass of 
each lepton or quark is less - and for the least massive 
particle in each family, much less - than the sum of 
the masses of the constituents. 

For simplicity the bound system would likely be 
comprised of two or three constituents. If the com- 
posite system were comprised of a spin-0 boson and 
a spin-1/2 fermion, then, in the language of non- 
relativistic physics, if all states have zero orbital an- 
gular momentum, all would have total angular mo- 
mentum or spin one half. For the relativistic equation 
discussed here, a similar mechanism allows strongly 
bound states only when the total angular momen- 
tum or spin is one half. If the composite system were 
comprised of three or more constituents, it is difficult 
to find a mechanism that would prevent higher-spin 
bound states. 

Many physicists in the 1970's studied the possibil- 
ity that quarks and leptons might be bound states of 
preons Q . Assuming the existence of only a few pre- 
ons, the existence of all quarks and leptons could be 
explained with each being a preon bound state. In 
the book The Trouble with Physics @, Smolin writes, 
"Unfortunately, there were major questions that the 
preon theories were not able to answer. These have to 
do with the unknown force that must bind the preons 
together into the particles that we observe. The chal- 
lenge was to keep the observed particles as small as 
they are while keeping them very light. Because preon 
theorists couldn't solve this problem, preon models 
were dead by 1980." 

Leptons interact gravitationally, weakly, and elec- 
tromagnetically. Among the three forces, the electro- 
magnetic interaction is the only one that might be able 
to provide the requisite strong binding. At sufficiently 
high energies, the electromagnetic interaction must, 
of course, be replaced by the electroweak interaction, 
but this added complication is not incorporated in the 
Bethe-Salpeter equation discussed here. The similari- 
ties of the charged lepton and quark mass spectra sug- 
gest that the same mechanism might be responsible for 
binding in all four families. Such a scenario could oc- 
cur if only one of the two quark constituents interacted 
strongly so that the pair would not. A heretofore un- 
known interaction might, of course, be responsible for 
the binding, but assuming the existence of such a force 
would represent a much more speculative approach. 

A few nonrelativistic and partially relativistic cal- 
culations suggest that electromagnetism might be able 
to create a bound state with the properties of quarks 
and leptons. Using the Schrodinger equation in two 
space dimensions, a charged quanta interacting with 
a charged magnetic dipole was shown to create bound 
states such that all low-energy states have zero or- 



bital angular momentum [lOj. The Dirac equation 
was used to show that strongly bound states can occur 
when a spin-1/2 quanta interacts with a charged mag- 
netic dipole |11| . and the Klein-Gordon equation in 
two space dimensions was used to show that a charged 
quanta interacting with a charged magnetic dipole can 
create bound states such that the energy gap increases 
between successively higher bound states [13] • 

But, as previously mentioned, by about 1980 most 
physicists had given up on the idea: experimen- 
tally there was no conclusive evidence that quarks 
or leptons were composite, a situation that has not 
changed significantly in the intervening thirty years, 
and theoretically no one could find a way to create 
bound states that both possessed the correct mass 
spectra and also were extremely small. The Bethe- 
Salpeter equation [l^ provided a theoretical frame- 
work for studying relativistic bound states. Also, by 
the late 1960's two, two-body, bound-state Bethe- 
Salpeter equations with unphysical interactions had 
been completely solved. However, in 1980 a gen- 
eral, numerical technique for solving the two-body, 
bound-state Bethe-Salpeter equation had not been de- 
veloped, and the high-speed computers required to 
solve the bound-state equation when the preons inter- 
act via minimal electrodynamics did not exist. With 
the development of super computers and a general 
method for solving the two-body, bound-state Bethe- 
Salpeter equation |14i] . the properties of relativistic 
bound states can now be determined, providing moti- 
vation to revisit preon models. 



II. POSSIBLE PREON MODEL OF LEPTONS 
AND QUARKS 

If the charged lepton family, the neutrino family, 
the negatively charged quark family, and the posi- 
tively charged quark family are each a bound state 
of a different combination of a spin-0 and a spin-1/2 
preon that interact via minimal electrodynamics, then 
all four families could be created from just two fermion 
preons and two boson preons. Each of the two boson 
preons can bind with each of the two fermion preons 
to create four different bound states or four families 
as indicated in Table 1. Also, as indicated in Table 
1, the sum of the charges of the two constituent pre- 
ons must equal the charge of the leptons or quarks in 
the family that the preons combine to create, yield- 
ing four equations for the charges of the four preons. 
Only three of the equations are independent. Possi- 
ble preon charges < 2 that are multiples of ±i are as 
follows: 
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TABLE I: Preon charges and constraints that the preon 
charges must obey 

Family Preon Preon Constraint 

Fermion Boson 
Charge Charge 



Electron gyi Qsi 

Neutrino 5/2 Qsi 

Negative Quarks g/i Qs2 

Positive Quarks g/2 Qs2 



qfi + Qsi = -1 
(7/2 + Qsi = 
(7/1 + Qs2 = -| 

9/2 + QS2 = I 



TABLE II: Possible charges of preon fermions and bosons 



Preon Fermion Charges 



Preon Boson Charges 



g/i 
1 
2 
3 
1 
3 

_ 4 
3 

_ 5 
3 

-2 



9/2 

2 

5 
3 
4 
3 

_ 1 
3 

_ 2 
3 



QS2 



III. INTRODUCTION TO THE 
BETHE-SALPETER EQUATION 

For relativistic, two-body, bound-state systems, a 
few analytical, Bethe-Salpeter solutions exist, but 
only in the limit that the binding energy equals the 
sum of the masses of the two constituent particles, 
implying the the energy of the bound state is zero. 
These solutions are not physical because they are cal- 
culated in the center-of-mass rest frame. But if a 
bound state has zero rest energy, it must travel at 
the speed of light. Therefore, zero-energy solutions 
are the zero-energy limit of physical solutions. In this 
limit the equation is always invariant under rotations 
in Minkowski space and so is separable. 

Analytical solutions have been obtained in the zero- 
energy limit, for example, for the Wick-Cutkosky 
Model [1^, [3| that consists of two scalars with ar- 
bitrary, nonzero masses bound by the exchange of a 
massless scalar. A variety of exact solutions have been 
obtained for two spinors with equal masses that are 
bound by the exchan ge o f a massless scalar [H, [l^ 
or a massless vector [23^. An exact solution has also 
been found for two spinors with unequal, nonzero 
masses bound by the exchange of a massless scalar 

Even when the energy is zero, most solutions are 
numerical. For example, such solutions have been ob- 
tained for two scalars with unequal, nonzero masses 
bound by the exchange of a massive scalar [2§| and for 
a scalar and a spinor with arbitrary, nonzero masses 



bound by the exch ang e of a massless scalar [23], |2J] or 
a massless vector [25[. When the energy is nonzero, 
all solutions are numerical. 

Only four, two-body, bound-state Bethe-Salpeter 
equations have been completely solved. That is, 
solved for arbitrary values of energy and arbitrary val- 
ues of the masses of the two constituents. (1) The 
Wick-Cutkosky Model consists of two scalars with ar- 
bitrary, nonzero masses bound by the exchange of a 
massless scalar [l^ M, Hi- (2) The Scalar-Scalar 
Model consists of two scalars with arbitrary, nonzero 
masses bound by the exchange of a massive scalar 
[13, SI m, [13 ■ (3) The Scalar Electrodynamics 
Model consists of a scalar and a fermion with arbi- 
trary, nonzero masses bound by the exchange of a 
massless scalar [l3|- (4) The Fermion- Scalar Model 
consists of a scalar and a fermion with arbitrary, 
nonzero masses bound by the exchange of a massive 
scalar f3l|. 

The interaction of the preons and the electromag- 
netic field is described by the renormalizable interac- 
tion Lagrangian, 



Lint = ■ -qfA'^ip^^-ip + iQs Af, 
+ QlA^A.cj,cb^ : . 



[0 {d^'(j)^)-{d^'^) (j)^] 
(2) 



The proposed preon model of leptons, quarks or both 
is a bound state of a spin-1/2 fermion with charge qf 
and mass m/ represented by the field tp ^J^d a scalar 
with charge Qs and mass mg represented by the field 
(f). The fermion and scalar fields interact minimally 
with the electromagnetic field A^. 

The Bethe-Salpeter equation will be used to study 
the properties of the bound state. Since every Feyn- 
man diagram contributes to the exact equation, some 
approximation must be made. Usually only the effect 
of the lowest-order diagram shown below is included. 
Because of the structure of the Bethe-Salpeter equa- 
tion, contributions from the diagram and iterations 
of the diagram, which form a ladder with the photon 
propagator forming the rungs, contribute. Thus the 
approximation is called the ladder approximation. 

In the ladder approximation the two-body, bound- 
state Bethe-Salpeter equation describing a spin-0 bo- 
son and a spin-1/2 fermion bound by minimal electro- 
dynamics is 

[(/ + eS)7o+p^7.-"^/] 



(27r)4 {p - qY -f le 

+ l\p^ + q^)}x{q)■ 



{7°bo 



qo 



2{i~l)E] 
(3) 



The equation has been written in the center-of-mass 
rest frame where the four-momentum of the bound 
state takes the form if = (i?, 0,0, 0). In the non- 
relativistic limit ^ = m//(m/ 4- "nis). 
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IV. STRONGLY BOUND SOLUTIONS 



9/ 



To solve the Bethe-Salpeter equation in the strong- 
binding (zero-energy) Umit, the singularity in the ker- 
nel of the Bethe-Salpeter equation is removed, and the 
equation is transformed from Minkowski to Euclidean 
space by making a Wick rotation ^] , which is always 
possible in the ladder approximation. 

In the zero-energy limit the angular dependence of 
the equation separates [32] when solutions are writ- 
ten as products of the functions and '^'"2^^ that 
are four-component hyperspherical harmonics in four- 
dimensional, Euclidean space-time and the functions 
F{\p\) and Gdpl) that depend only on the magnitude 
of the Euclidean four-momentum \p\, 



FIG. 1: Lowest-order Feynman diagram for a spin-1/2 field 
ijj with charge Qf and a scalar field (j) with a charge Qs 
interacting via minimal electrodynamics with the photon 
field A^". 



If the binding were non-relativistic, the properties 
of the bound states could be calculated using the 
Schrodinger equation. 



X(*p°,p) = F{\p\)¥+\cos e^,92, 

+ G(|p|)*^+)(C0S 01,02, 



(4) 



There are a number of notable differences between 
the bound-state Schrodinger and Bethe-Salpeter equa- 
tions: a) Because the energy E appears multiple times 
in the Bethe-Salpeter equation ([3]), a Hamiltonian 
does not exist for the relativistic problem, b) Because 
E appears multiple times in the Bethe-Salpeter equa- 
tion, the equation is solved by specifying the energy 
and solving for the coupling constant q/Qs/ (47r) as an 
eigenvalue. The process for solving the Schrodinger 
equation ^ is the reverse: the coupling constant is 
specified, and the energy E is solved for as an eigen- 
value. The two procedures yield equivalent informa- 
tion, (c) There is no action at a distance for the 
Bethe-Salpeter equation since the interaction is co- 
variant, (d) The Bethe-Salpeter is an integral equa- 
tion; therefore, boundary conditions are incorporated 
in the equation, (e) The Bethe-Salpeter equation is 
separable when energy E ^ and is almost always 
nonseparable when E ^ 0. 

The author has developed a systematic method [ij] 
for solving the two-body, bound-state Bethe-Salpeter 
equation for arbitrary values of energy. However, to 
show how specific boundary conditions allow for the 
possibility that strongly bound solutions exist when 
the coupling constant g/(3s/(47r) is on the order of the 
fine structure constant a, here attention is restricted 
to the strong binding limit E = 0. 



(5) 



(+) 



The four-component hyperspherical harmonics ^'j^ 

and 5* 2^'' in ([5]) are analogous to the two-component 

spherical harmonics [33] that are used to sepa- 
rate the Dirac equation with a spherically symmet- 
ric potential. Using hyperspherical harmonics in four- 
dimensional space as part of the basis functions both 
decreases the number of integrations that must be per- 
formed numerically and prevents the appearance of a 
logarithmic singularity in the kernel. The separated, 
zero-energy Bethe-Salpeter equation is 



[IpP + (1 + A)2] [(1 - A)F(+)(|p|) + \p\G(+H\p\) 



(2^)4 
QfQs 

(2^)4 



9l'Ai'Vi/2(bU'?l)G^+'(l9l)d|'Zl 







q\%%/,{\p\,\q\)G^^H\qMq\, (6a) 



[\p\^ + (1 + A)2] [-|p|F(+)(b|) + {1- A)G(+)(|p|) 



(2^)4 1^1 



q\'AZi,/,{\pl\q\)F^+\\q\)d\q\ 



(27r)-^o 



9K'+i/2(bU'?l)^^+Hkl)d|g|. (6b) 



In the above equation the index fci = l/2,3/2,... and 

' 111 yi+l 



^^'H\p\,\q\) 



|f&%^^ if kl<H 



, I IPI \n + i 



(7) 



Dimensionless variables have been introduced in (6) 
by defining m/ = m(l — A),ms = to(1-|-A),p' = p/m 
and q' = q/m. Then primes have been omitted for all 
momentum variables since all are dimensionless. 



Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



5 



Eq. (6) is solved by expanding the functions i^(|p|) 
and G(|p|) in terms of a finite set of basis functions 
that (very nearly) obey the boundary conditions that 
F{\p\) and G{\p\), respectively, obey. 

basis function = convergence function x spline (8) 

The convergence function depends on the magnitude 
of the four- momentum \p\ and obeys the boundary 
conditions that the solution obeys. A general method 
for analytically calculating boundary conditions for 
Bethe-Salpeter equations has been developed 
Typically, as is the case for the Bethe-Salpeter equa- 
tion being discussed here, the behavior of the solution 
at small momenta |p| can be determined exactly; how- 
ever, at large |p| boundary conditions can sometimes 
only be determined within ranges, complicating the 
calculation. If the solution obeys the boundary 

conditions 



#1 #2 



F{\p\) bK" 



n\p\) 



\p\- 



\p\- 



(9) 



then the convergence function for is given by 



convergence function 



bl 



/o 



(IpI + constant)/o+/°< 



(10) 



Cubic splines _Bj(|p|) 34] are defined on four con- 
tiguous intervals, and the value of \p\ at the bound- 
ary of each interval is called a "knot" that is denoted 
by Tj. Within each interval a cubic spline is a cubic 
polynomial with coefficients chosen so that the cubic 
spline and its first two derivatives are continuous, im- 
plying that a cubic spline and its first two derivatives 
vanish at the beginning of the first interval and at 
the end of the fourth. The single, remaining, unspec- 
ified coefficient is fixed by an overall normalization 
condition. The spline Bi{\p\) begins at the knot Ti 
and terminates at T'i+4. A major advantage of using 
splines, instead of functions defined on the entire in- 
terval < IpI < 00, is that, by concentrating knots in 
the region where the solution is changing most rapidly, 
more splines are available where they are most needed 
to represent a solution. Additionally, it is much faster 
to integrate numerically over splines because they are 
nonzero only on a finite interval. 

Because the convergence function obeys the bound- 
ary conditions exactly at large and small \p\, and be- 
cause the splines also depend on \p\, the basis func- 
tions very nearly obey the boundary conditions, but 
they do not obey the boundary conditions exactly. 
However, linear combinations of these basis functions 
obey the boundary conditions exactly and almost al- 
ways yield convergent series expansions for solutions. 

The zero-energy equation is discretized and then 
solved by converting it into a generalized matrix eigen- 
value equation usin g th e Rayleigh-Ritz-Galerkin vari- 
ational method ISa, Sol • 




FIG. 2: Cubic spline i3i(|p|) plotted as a function of the 
magnitude of the four-momentum \p\. 



Typically strongly bound solutions of the Bethe- 
Salpeter equation exist only when the coupling con- 
stant qfQs/{4:'K) is on the order of or greater than 
unity. If solutions exist for g/(5s/(47r) ^ a, then at 
least one integral in (6) must be very large. Two sets 
of boundary conditions [2^ could yield one or more 
large integrals, and in both cases the integrals could 
be large when the coupling constant is small only if the 
spin of the bound state is one half. If leptons, quarks 
or both are composite, this mechanism would explain 
why leptons, quarks or both all have spin one half. 
Here attention is restricted to the following boundary 
conditions: 



F{\P\) 
Gi\p\) 



In (11), fci 



\p\^oo 



IpI- 



Foo \P\ 
Goo \P\ 



'(fcl+i) 



-(fcl + l+eg) 



J, J 



1,. 



and Foo and Gn 



(11a) 
(lib) 
are 



constants. The parameter eg satisfies the conditions 
< eg < 2 if fci = 1/2, and < eg < 2 if fci > 3/2. 

For solutions satisfying the boundary conditions 
(11), in the limit of large \p\, the Betlic-Salpcter equa- 
tion (|6ap decouples from (j6bp and takes the form 



Goo\p\ 



2{ki + 1) 



^ qjQs 1 

■\p\ 



(fci + 5)(fci + |) 

|p|*:i+3/2 



bl 



-(fci+3/2) 



d|g||g 



fei+7/2 



G{\q\) 







^1 

|p|fei-l/2 
fcl + ^ 



d|g|M 



Go 



d|g||g 



-(2fei+eg) 



(12) 



Only the final integral in can be large and then 
iff fci = 1/2 and eg <^ 1. Since fci = j,j + I, . . . , 
and fci = 1/2 for all strongly bound states with small 
coupling constants, all such states have total angular 
momentum or spin one half. 
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Using expressions for the integrals in (fT^ that are 
valid in the limit of large |p| [2J], the following ex- 
pression is obtained for the coupling constant when 
ki = 1/2: 

qjQs _ 7r eg(2-eg)(2 + eg) 

For small values of eg, the coupling constant 
g/(5s/(47r) is small. 

It is particularly difficult to obtain solutions when 
a fermion and scalar interact via minimal electrody- 
namics instead of through the exchange of a massless 
or massive scalar. As can be seen from the depen- 
dence of the coupling constant on Cg, the solutions go 
to zero at large momenta sufficiently slowly that the 
behavior of solutions at large momenta has a very 
significant effect on the coupling constants. 



V. SUMMARY 

Three major points are made in the talk: a) It 
might be possible to describe leptons, quarks or both 
as highly relativistic bound states of a spin-0 and spin- 
1/2 constituent bound by minimal electrodynamics, 
b) Strongly bound states with coupling constants on 
the order of the electromagnetic fine structure con- 
stant a are allowed by the boundary conditions and 
likely exist, c) All such strongly bound states would 
have spin one half, as do the leptons and quarks. 
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